LAMMPS (18 Feb 2013)
# Pour 2d granular particles into container

dimension	2
atom_style	sphere
boundary	f fm p
newton		off
communicate	single vel yes

region		reg block 0 100 0 50 -0.5 0.5 units box
create_box	1 reg
Created orthogonal box = (0 0 -0.5) to (100 50 0.5)
  4 by 1 by 1 MPI processor grid

neighbor	0.2 bin
neigh_modify	delay 0

# IMPORTANT NOTE: these values of stiffness (4000) and timestep (0.001)
# are used in this example file to produce a quick simulation and movie.
# More appropriate values for realistic simulations would be
# k = 4.0e5 and dt = 0.0001, as in bench/in.chute (if it were Hertzian).

pair_style      gran/hertz/history 4000.0 NULL 100.0 NULL 0.5 0
pair_coeff	* *

timestep	0.001

fix             1 all nve/sphere
fix		2 all gravity 1.0 spherical 0.0 -180.0
fix		xwalls all wall/gran 4000.0 NULL 100.0 NULL 0 1 		xplane 0 100
fix		ywalls all wall/gran 4000.0 NULL 100.0 NULL 0 1 		yplane 0 NULL

region          slab block 1.0 99.0 30 34.5 -0.5 0.5 units box
fix             ins all pour 1000 1 4767548 vol 0.4 10 	        diam 0.5 1.0 region slab
Particle insertion: 224 every 3000 steps, 1000 by step 12001

fix             3 all enforce2d

compute		1 all erotate/sphere
thermo_style	custom step atoms ke c_1 vol
thermo		1000
thermo_modify	lost ignore norm no
compute_modify	thermo_temp dynamic yes

#dump		id all atom 250 dump.pour

#dump		1 all image 500 image.*.jpg type type #		zoom 1.6 adiam 1.5
#dump_modify	1 pad 5

run		25000
Memory usage per processor = 9.54554 Mbytes
Step Atoms KinEng 1 Volume 
       0        0           -0            0         5000 
    1000      224    185.87695            0         5000 
    2000      224    348.08417            0         5000 
    3000      224     565.2162            0         5000 
    4000      448    1042.9716            0         5000 
    5000      448    1542.0614            0         5000 
    6000      448    1766.5136    41.192729         5000 
    7000      672    1846.4596    67.352592         5000 
    8000      672    1814.3429    92.898167         5000 
    9000      672    1662.3643    111.88767         5000 
   10000      896    1652.8654    91.325153         5000 
   11000      896    1541.4417    70.862864         5000 
   12000      896    1474.3165    61.956769         5000 
   13000     1000    1463.9007    59.436066         5000 
   14000     1000     1355.973    51.999639         5000 
   15000     1000    1130.5891     49.13995         5000 
   16000     1000    787.89487     46.05549         5000 
   17000     1000    462.83858    40.678165         5000 
   18000     1000    325.64118    32.743319         5000 
   19000     1000    155.22607    30.549487         5000 
   20000     1000    43.494925    22.262581         5000 
   21000     1000    23.694914    15.061018         5000 
   22000     1000    20.580126    10.477019         5000 
   23000     1000     16.06558    7.3034265         5000 
   24000     1000    12.026429    6.0665971         5000 
   25000     1000    11.818403    5.4542268         5000 
Loop time of 0.908638 on 4 procs for 25000 steps with 1000 atoms

Pair  time (%) = 0.304472 (33.5086)
Neigh time (%) = 0.0796848 (8.76969)
Comm  time (%) = 0.0719426 (7.91763)
Outpt time (%) = 0.000483155 (0.0531735)
Other time (%) = 0.452056 (49.7509)

Nlocal:    250 ave 252 max 248 min
Histogram: 1 0 1 0 0 0 0 1 0 1
Nghost:    18 ave 25 max 12 min
Histogram: 1 1 0 0 0 0 1 0 0 1
Neighs:    604.25 ave 614 max 589 min
Histogram: 1 0 0 0 1 0 0 0 0 2

Total # of neighbors = 2417
Ave neighs/atom = 2.417
Neighbor list builds = 2023
Dangerous builds = 0
